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Abstract 

Background: Increasing structural and biochemical evidence suggests that post-translational methionine oxidation 
of proteins is not just a result of cellular damage but may provide the cell with information on the cellular oxidative 
status. In addition, oxidation of methionine residues in key regulatory proteins, such as calmodulin, does influence 
cellular homeostasis. Previous findings also indicate that oxidation of methionine residues in signaling molecules 
may have a role in stress responses since these specific structural modifications can in turn change biological 
activities of proteins. 

Findings: Here we use tandem mass spectrometry-based proteomics to show that treatment of Arabidopsis thaliana 
cells with a non-oxidative signaling molecule, the cell-permeant second messenger analogue, 8-bromo-3,5-cyclic 
guanosine monophosphate (8-Br-cGMP), results in a time-dependent increase in the content of oxidised methionine 
residues. Interestingly, the group of proteins affected by cGMP-dependent methionine oxidation is functionally 
enriched for stress response proteins. Furthermore, we also noted distinct signatures in the frequency of amino acids 
flanking oxidised and un-oxidised methionine residues on both the C- and N-terminus. 

Conclusions: Given both a structural and functional bias in methionine oxidation events in response to a signaling 
molecule, we propose that these are indicative of a specific role of such post-translational modifications in the direct or 
indirect regulation of cellular responses. The mechanisms that determine the specificity of the modifications remain to 
be elucidated. 

Keywords: Methionine oxidation, Reactive oxygen species, 3,5-cyclic guanosine monophosphate, 
Tandem mass spectrometry-based proteomics, Arabidopsis thaliana 



Findings 

The debate of whether methionine (Met) oxidation of 
proteins is a purely chemical consequence of cellular 
oxidative damage or a protective mechanism against oxi- 
dative damage, or indeed a post-translational modifica- 
tion that can act as a specific cellular signal and/or 
response, is ongoing [1-3]. To shed light on this question 
we treated Arabidopsis suspension culture cells with the 
cell permeant second messenger analogue 8-bromo 3,5- 
cyclic guanosine monophosphate (8-Br-cGMP). Cyclic 
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GMP has a signaling role in many plant responses, in- 
cluding responses to light [4], hormones [5-8], signaling 
peptides [9], salt and drought stress [10,11], ozone, and 
defence responses [12-14]. Given that cGMP is not an 
oxidising agent and does not induce protein Met oxidation 
in vitro (Additional file 1), we tested if cGMP causes pro- 
tein oxidation in vivo. To this end we used an OxiSelect™ 
Intracellular ROS Assay Kit (Cell Biolabs, Inc.) and show 
that cGMP can cause protein oxidation (Figure 1). We 
therefore conclude that any protein Met oxidation event 
resulting from cGMP treatment is most likely the result of 
direct or indirect cellular processes. To further character- 
ise cGMP-dependent Met oxidation in vivo, a proteomic 
analysis was performed on A. thaliana (ecotype Col-0) cell 



© 2013 Marondedze et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the 
Creative Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly cited. 



Marondedze et al. Cell Communication and Signaling 2013, 11:1 
http://www.biosignaling.eom/content/1 1 /I /I 



Page 2 of 7 



20000 
18000 
16000 
14000 
12000 
10000 
8000 
6000 
4000 
2000 
0 



Control 

30' 
60' 



, il ''J. 
I I I I 



OmMcGMP IOhMcGMP 50hMcGMP 10mMH 2 O 2 

Figure 1 Protein oxidation assay. OxiSelect™ Intracellular ROS 
assay kit (Cell Biolabs, Inc. San Diego, CA) was used in the in vivo 
oxidation experiments according to the assay protocol provided by 
the manufacturer. Cultured Arabidopsis (Col-0) cells were placed in a 
black bottom 96-well cell culture plate for 2 h in a shaking 
incubator. The 2 / ,7 / -dichlorofluorescein diacetate/media solution was 
added to the cells prior to incubation at 37°C for 1 h. The 
dye-loaded cells were then treated with 10 uM or 50 uM of cGMP or 
H 2 0 2 . Fluorescence in the cells was measured at 30 and 60 min 
post-treatment at 480/530 nm using a PHERAstar FS microplate 
reader (BMG Labtech GmbH, Germany) and the values plotted. Each 
bar represents data from 3 biological replicates (n = 3), the bars are 
the standard errors. Treatment with 8-Br-cGMP at the final 
concentration of 50 uM induces statistically significant differences of 
the means at p = 0.05 using a two-sample f-test. 



suspension culture grown in Murashige and Skoog medium 
[15] following treatment with 8-Br-cGMP at the final con- 
centration of 10 \iM. Three biological replicate samples 
were collected at 0, 30 and 60 minutes post-treatment 
Total soluble proteins were extracted [16] and processed 
for tandem mass spectrometric identification of peptides 
containing oxidised Met residues (for methods see legend 
to Figure 2). Since our experimental protocol included a 
Ti0 2 enrichment step, usually applied for enrichment of 
phosphopeptides, we also tested to what extent peptides 
containing oxidized Met residues are enriched in samples 
subjected to the Ti0 2 -based enrichment step either in the 
presence of absence of DHB (2,5-dihydroxybenzoic acid). 
In both cases the enrichment led to significant increase in 
ratio of spectra assigned to oxidised Met peptides to all 
assigned spectra (Table 1). Moreover, the presence of 
DHB in the Ti0 2 -based enrichment step enhanced further 
increase in the number of oxidized peptides identified as 
compared to oxidized Met peptides enriched in the ab- 
sence of DHB. This is consistent with a report that shows 
that Met oxidised peptides co-enrich with phosphopep- 
tides because the affinity for the Ti0 2 (in the presence of 
DHB) is stronger in oxidised as compared to non-oxidised 
isoforms [17]. 



In our proteomic analysis we considered a peptide as 
containing oxidised Met residue when it was identified 
with high confidence (> 95%) in at least two biological 
replicates. A total of 385 cGMP-dependent methionine 
oxidised proteins were identified (Additional file 2, tab 
"AF1"). Assigned spectral counts (Additional file 2, tab 
"AF2") were used to estimate the relative ratio of peptides 
containing oxidized Met residue(s) as compared to total 
number of peptides identified in the sample. 

An example of a tandem mass spectrometry result 
demonstrating oxidative modification of Ti0 2 -enriched 
peptides extracted from 8-Br-cGMP-treated cells is 
shown in Figure 2. Peptides containing single oxidised 
Met residue show an increase in mass to charge (m/z) 
ratio of 15.9994 that corresponds to the average mass of 
an oxygen atom. For example, the peptide fragment 
(DHDKPIQQVIAEMTDGGVDR) of the alcohol de- 
hydrogenase 1 (AT1G77120) in non-oxidised form has 
the m/z ratio of 850.3723 (Figure 2A), while after oxida- 
tion of Met residue, the m/z ratio shifts to 866.3673 
(Figure 2B). 

Further, we identified peptides with oxidised Met that 
occurred in all three biological replicates at different 
time points. We noted an increase in the total number 
of peptides containing residues of oxidised Met after 
cGMP treatment from 221 to 633 and then 1451 at 0, 30 
and 60 minutes, respectively (Figure 3A and Additional 
file 2, tab "AF2"). These numbers represent 1.4%, 19.4% 
and 13%, respectively, of the total number of peptides 
identified at each time point. Thus, the percentage of 
Met oxidised peptides identified is the highest at 30 min- 
utes. In addition, the numbers of oxidised Met peptides 
detected at each time-point suggest that the total num- 
ber of oxidised Met residues increased nearly 3-fold dur- 
ing the first 30 minutes of treatment and 7-fold after 
60 minutes of treatment (Additional file 2, tab "AF2"). 
Of these redundant peptide fragments containing oxi- 
dized Met, 14 at 0 minutes, 113 at 30 minutes and 288 
at 60 minutes were unique for each time-point. The total 
Met oxidised peptides correspond to 34 (at 0'), 136 
(at 30') and 281 (at 60') Arabidopsis proteins, from 
which 10, 94, and 224 identified oxidized Met proteins 
are unique for each time-point, respectively (Figure 3B). 
This finding implies either that cGMP-dependent Met 
modifications are reversible and/or that some of the 
modified proteins have undergone proteolysis. 

Gene ontology (GO) analysis (Fatigo + ; http://babelomics3. 
bioinfo.cipf.es/) [18] of the unique proteins containing oxi- 
dised Met residues was undertaken. The result shows a 
significant enrichment of GO terms (adjusted p-value < 
l.OOe" 02 ) in categories including 'response to stress', 'response 
to abiotic stimuli', 'response to oxidative stress', and 'response 
to oxygen and ROS metabolic process' (see Additional file 2, 
tab "AF3"). Moreover, the gene number in these GO 
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Figure 2 MS/MS spectra of ADH1 containing non-oxidised (A) and oxidised (B) methionine residues. Three biological replicates of 10 uM 
8-Br-cGMP-treated cells and H 2 0 mock treated controls were collected at 0, 30 and 60 min. Proteins were precipitated using 10% (w/v) 
trichloroacetic acid in acetone, re-solubilised in 7 M urea, 2 M thiourea and 4% (w/v) CHAPS, reduced, alkylated and trypsin digested. Peptides 
were fractionated by cation exchange chromatography. Methionine oxidised peptides were enriched using Ti0 2 beads and re-suspended in 5% 
(v/v) acetonitrile and 0.1% (v/v) formic acid prior to identification and quantitation by LTQ Orbitrap coupled with a nanoelectrospray ion source. 
Peptides (5 uL) were injected onto a 50 mm x 0.3 mm Magic C18AQ column. The top 10 precursor ions were selected with a resolution of 60,000 
for fragmentation using normalized collision-induced dissociation set at 35.0. Spectra were searched against TAIR using MASCOT, with a precursor 
mass tolerance of 10 ppm, a fragment ion mass tolerance of 0.3 Da, one missed cleavage, carbamidomethyl cysteine residues as fixed 
modification and oxidation and dioxidation of methionine residues as variable modifications. Proteins with a score > 95% were considered 
positively identified (corresponding score < 31). Spectra were further processed with the Scaffold™ software using the "Trans-Proteomic Pipeline" 
algorithm (threshold 95%). Oxidised Met residues showed an increase in mass/charge ratio (m/z) of 15.9994. Arrows show Met residues at 
position 13 in the fragment DHDKPIQQVIAEMTDGGVDR of AT1G77120 before oxidation (m/z ratio 850.3723) (A) and after oxidation (m/z ratio 
866.3673) (B). 



categories also increased over time (Figure 4). It therefore 
appears that the application of a membrane-permeable 
analogue of the second messenger cGMP induces Met oxi- 
dation in a set of proteins that are over-represented in spe- 
cific functional groups. 

Proteins enriched in the GO categories: 'response to oxi- 
dative stress', 'response to ROS', 'response to oxygen', and 
'ROS metabolic process' showed four main patterns. 
(1) Loss of peptides or proteins containing oxidised Met 
residue(s) and reduction in the number of peptide copies 
after treatment, e.g. glyceraldehyde-3-phosphate dehydro- 
genase (AT1G13340; AT3G04120) and 60S acidic ribo- 



somal protein (AT2G27720). Loss of Met oxidised peptides 
can occur due to degradation of copies of proteins with 
oxidised Met residues or the reduction of oxidised Met 
residues. (2) Increase in the number of oxidised peptide 
fragments after treatment, e.g. peroxidase (AT2G22420; 
AT4G08770). This may imply that cGMP either indirectly 
induces over- expression of specific proteins with peptide 
fragments susceptible to oxidation or preferentially 
induces oxidation of Met residues in specific proteins 
without necessarily inducing transcription and/or transla- 
tion. (3) New peptides detected with Met residues oxi- 
dised after treatment, e.g. ATP synthase (AT5G08670), 



Table 1 Enrichment of methionine oxidized peptides (oxMet) using Ti0 2 with and without DHB 

Without enrichment Ti0 2 enrichment without DHB Ti0 2 enrichment with DHB 

R1 R2 R3 R4 R1 R2 R3 R4 R1 R2 R3 R4 

Assigned spectra 123109 67688 42897 89362 38403 56545 14450 4288 28092 11880 8567 3812 

Assigned spectra oxMet pep. 3073 288 49 3930 15803 4273 6303 424 26913 4370 6586 1318 
%oxMet 2.5 0.4 0.1 4.4 41.2 7.6 43.6 9.9 95.8 36.8 76.9 34.6 



+ The cell suspension culture was treated with 10 uM H 2 0 2 for 30 minutes. The experimental procedures are as described in the legend to Figure 2. 
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Figure 3 Quantitative representation of peptides (A) and 
proteins (B) containing oxidised Met residues identified by 
mass spectrometry. The total numbers of peptides and proteins 
containing at least one oxidised Met residue identified by LC-MS/MS 
in protein samples extracted from A tholiono cells treated with 
10 uM 8-Br-cGMP and collected at 0, 30 and 60 min post-treatment 
were analysed. Proteins matching the peptides were identified by 
searching against TAIR10 database using MASCOT and TPP 
algorithm, and only proteins within a 95% confidence threshold 
were considered. The total number of peptides (A) or proteins 
(B) are indicated outside the Venn diagram (in bold), the numbers 
inside are the unique peptides or proteins identified at each time 
point. Fourteen, 113 and 288 peptide fragments containing at least 
one oxidised Met appear only at either 0, 30 or 60 min. This 
corresponds to an increase in the total number of time point 
specific Met oxidised proteins from 10 to 94 and then to 224 (B). 



and (4) multiple Met residues becoming oxidised after 
treatment, e.g. in the heat shock protein 70 (AT3G 12580) 
(Additional File 2, tab "AF1"). 

A total of 299 proteins (78% of all unique proteins con- 
taining oxidised Met residues) comprised only one oxidised 
Met residue, while 86 proteins (22% of all unique proteins 
containing oxidised Met residues) had>2 oxidised Met 
residues (Additional file 2, tab "AF4"). Of those multi- 
oxidised 86 proteins 10 proteins had > 4 Met residues oxi- 
dised (Additional file 2, tab "AF5")- The low expression of 
osmotically responsive gene V (AT1G56070), 'tubulin 9' 
(AT4G20890), 'heat shock protein 70' (AT5G02500) and 
male gametophyte defective V (AT2G21870) contain the 



largest number of oxidised residues, representing 16%, 
23%, 29% and 50%, respectively, of their total Met content 
(Table 2 and Additional file 2, tab "AF5"). It is of interest 
that heat shock proteins AT5G02500, AT3G 12580 and 
AT5G09590 each had at least four oxidised Met residues. 

The protein with the highest number of Met oxidised 
peptide fragments in response to cGMP is the methyles- 
terase PCR A (AT1G 11580) with 235 oxidised fragments 
representing ~ 30% of the identified peptides (Table 3 and 
Additional file 2, tab "AF5"). Among the proteins with 
peptides show the greatest extend of Met oxidation, per- 
oxidase superfamily protein (AT2G22420) showed the 
highest ratio of Met oxidised peptides to all identified pep- 
tide fragments of the protein (114 out of 160), represent- 
ing « 70% of total peptides detected and assigned to this 
protein. Localised in the cytosol, this family of proteins is 
mainly involved in oxidation-reduction processes and 
responses to oxidative stress. In addition to ROS being 
generated by activation of a plasma membrane NADPH 
oxidase [19], an extracellular cell wall peroxidase is also 
involved in the biosynthesis of H 2 0 2 [20] and plays an im- 
portant role in plant resistance to pathogens [21]. A previ- 
ous study using antisense expression of a French bean 
peroxidase cDNA in Arabidopsis showed a reduction in 
mRNA level of peroxidase (AT3G49120), which in turn 
led to a reduction in the oxidative burst and eventually a 
reduced resistance to fungal and bacterial pathogens [22]. 
The authors hypothesized that peroxidases have a role in 
sustaining and/or initiating ROS that signal early defence 
responses in plants. Given that peroxidases themselves are 
highly susceptible to Met oxidation, and that this modifi- 
cation limits their activity, it suggests a negative feedback 
mechanism and links Met oxidation to the control of cel- 
lular redox balance. It has also been reported that 0 3 and 
NO, both oxidising agents, induce transcriptional activa- 
tion of scavenger-encoding genes, like alternative oxidase 
and glutathione peroxidase [12]. Furthermore, oxidation of 
Met residues has been shown to target both specific func- 
tional domains and consequently modify functional char- 
acteristics, e.g. in cytochrome c [23] and peroxidases [24], 
as well as regions outside functional domains, like in the 
nascent polypeptide-associated complex where modifica- 
tions do not seem to affect functionality. 

Given that Met oxidation can profoundly alter cellular 
responses, the question is if there is evidence for site se- 
lectivity and if so, what determines it. In order to ad- 
dress this question, we have subjected all 385 proteins 
with one or several oxidised Met to further analysis and 
noted the following. Firstly, the average Met frequency 
in the set of proteins containing oxidised Met is 0.027 as 
compared to the other amino acids (AA) in the complete 
proteome where it is 0.025 [25], and hence is nearly the 
same; we also noted that none of the 50 Arabidopsis 
proteins with the highest frequency of Met occurrence 
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Figure 4 Enriched gene ontology (GO) terms (adjusted p-value< I.OOe" 02 ). Unique proteins containing oxidised Met identified in 
cGMP-treated samples were used to search for GO term enrichment using Fatigo + (http://babelomics3.bioinfo.cipf.es/) and selected significantly 
enriched terms are represented. A significant increase in the enrichment of GO terms over time was observed as well as an increase in the 
number of genes in these GO categories. 



(0.09) contained any oxidised Met residues. This indi- 
cates that Met frequencies in proteins per se are not a fac- 
tor that determines increased oxidation. Secondly, of the 
385 proteins containing at least one modified Met residue, 
only eight residues were on the N-terminus, and of the 
150 double Met (-MM-) residues, five proteins had both 
residues oxidised and eight proteins - only one. Thirdly, of 
the 575 oxidised Met residues detected, 75 had a glutamic 



acid (Glu) and 68 had an aspartic acid (Asp) as an imme- 
diate C-terminal neighbour and this bias is likely, at least 
in part, due to preferential enrichment of these AAs by 
Ti02 [26], however we also find the uncharged alanine 
(Ala: 47) enriched on the C-terminus. The most frequent 
N-terminal neighbours are glutamic acid (Glu: 67) and 
again the uncharged alanine (Ala: 65). The least frequent 
C- and N-terminal neighbours are tryptophan (Trp: 2; 0, 



Table 2 cGMP-dependent proteins with multiple methionine residues oxidised and their enriched GO terms 



Accession No. 


Description 


*OxMet/total 
Met residues 


Annotations 


Enriched 
GO term(s) 


AT1G56070 


Low expression of osmotically responsive gene 1 


5/38 


RC, translational elongation 


RTS 


AT2G21870 


Male gametophyte defective 1 


5/10 


Co 2+ , Cu 2+ and Zn 2+ binding 




AT4G20890 


Tubulin 9 


5/22 


GTP binding, GTP catabolic process, RCa 




AT5G02500 


Heat shock protein 70 


5/17 


Responses to heat, bacterium, H 2 0 2 and HLI, 
RCa, RTS, RV 


RTS, RV 


AT3G 12580 


Heat shock protein 70 


4/16 


Responses to heat, bacteria, H 2 0 2 , HLI, RCa, 
RTS, RV 


RTS, RV 


AT3G27240 


Cytochrome CI family 


4/8 


Heme binding, Fe 2+ binding 




AT4G02930 


GTP binding elongation factor Tu family protein 


4/14 


ATP, Co 2+ and Zn 2+ binding, RCa 




AT1G24310 


Unknown protein 


4/12 






AT1G27390 


Translocase outer membrane 20 


4/8 


Metal ion binding 




AT4G 10480 


Nascent polypeptide-associated complex 


4/10 


Transcription regulation and mitochondrial 
translocation 





HLI- high light intensity, RC- response to cold, RCa- response to cadmium ion, RTS- response to temperature stimulus (GO:0009266), RV- response to virus 
(GO:0009615), adjusted p < 1.78e" 2 . *Ox- oxidised. 
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Table 3 Highly methionine oxidised proteins after cGMP treatment 



Accession 
No. 


Description 


% of oxidised 
fragments 


Annotations 


Enriched GO 
term(s) 


AT2G22420 


Peroxidase superfamily 
protein 


71.3 


Oxidation-reduction process, response to oxidative stress 




AT5G 12250 


(3-6-tubulin 


62.6 


Response to salt stress, RC 




AT4G 10480 


NAC 


50.9 


Transcription regulation and mitochondrial translocation 




AT3G 15950 


NAI2 


41.8 


Response to salt stress 




AT2G05710 


Aconitase 3 


35.0 


Response to ABA stimulus, oxidative and salt stress, RCa 




AT5G02500 


Heat shock protein 70 


33.6 


Defense response to bacteria and fungus, response to heat, 
and virus, RCa, RC 


RAS, RS 


AT2G21660 


Cold, circadian rhythm and 
RNA binding 2 


32.7 


Regulation of stomatal movement, response to osmotic and salt 
stress, RCa, RC 


RAS, RS 


AT1G1 1580 


Methylesterase PCR A 


30.8 


Metabolic process, negative regulation of catalytic activity 




AT1G77120 


Alcohol dehydrogenase 


18.7 


Cellular respiration, oxidation-reduction process, response to hypoxia, 


RAS, RS 



osmotic stress and salt stress, RCa 



AT1G53240 Mitochondrial malate 16.9 Oxidation-reduction process, defense response to bacteria, response 

dehydrogenase to salt stress, RCa, RC 

NAC Nascent polypeptide-associated complex, NAI2- similar to TSK-associating protein 1, RC Response to cold, RCa Response to cadmium ion, RAS Response to 
abiotic stimulus (GO:0009628), RS Response to stress (GO:0006950), adjusted p <4.15e~ 2 . 



respectively) and cysteine (Cys: 1; 2, respectively). We fur- 
ther compared the observed frequency of amino acids 
flanking Met positions to their theoretically expected 
values implied from the analysis of overall AA frequencies 
in the entire Arabidopsis proteome [25]. This process 
involved counting all of the occurrences of AAs in the 
proteome and establishing their relative frequency in the 
proteome, and the analysis was done in Matlab (Version 
R2010b). We note that, while under-represented (i.e. the 
observed frequency in positions flanking Met is lower than 
their average in the entire proteome), leucine and serine 
are the most frequent flanking AAs. In turn, Met flanking 
a Met is 35% over-represented on the N-terminal side and 
25% over-represented on the C-terminal side. In contrast, 
oxidised Mets have much reduced relative preference for 
Met (-14% on N- terminus and -44% on C-terminus), albeit 
based on a very limited sample. Other under-represented 
flanking AAs are cysteine (Cys) and proline (Pro). 

The rapidly evolving field of redox proteomics provides 
new evidence supporting the notion that oxidation of Met 
residues may have a great impact on protein activity, regu- 
lation of biochemical pathways and cellular function in 
response to changing environmental conditions. This is 
consistent with the observed Met oxidation accumulation 
in plants under low temperature conditions and the fact 
that plant methionine sulfoxide reductase (MSR) confers 
increased tolerance to freezing [27]. It is also conceivable 
that the differential oxidation footprint of Met is a re- 
sult of different susceptibility depending e.g. on the con- 
formation of the protein or on differential access for repair 
of the MSR to different proteins or protein domains. 
Furthermore, our results are an indication that many of the 
Arabidopsis proteins involved in modulating the level of re- 
active oxygen species (ROS), including ROS-scavenging 



and ROS-producing proteins [19], may - at least in part - 
be regulated by oxidation of their Met residues. 

Additional files 



Additional file 1: The file contains the list of methionine oxidised 
peptides and their assigned spectral counts following in vitro cGMP 
treatment experiment. 

Additional file 2: The file contains the lists of all methionine 
oxidised proteins and fragments and the result of the gene 
ontology analysis. 
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